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We study excitonic states of an atomic impurity in a Fermi gas, i.e., bound states consisting of the impurity 
and a hole. Previous studies considered bound states of the impurity with particles from the Fermi sea where 
the holes only formed part of the particle-hole dressing. Within a two-channel model, we find that, for a wide 
range of parameters, excitonic states are not ground but metastable states. We further calculate the decay rates 
of the excitonic states to polaronic and dimeronic states and find they are long lived, scaling as oc 
and oc (Aft))^. We also find that a new continuum of exciton-particle states should be considered alongside 
the previously known dimeron-hole continuum in spectroscopic measurements. Excitons must therefore be 
considered as a new ingredient in the study of metastable physics currently being explored experimentally. 


I. INTRODUCTION 

An impurity interacting strongly with a background gas is 
a paradigm of many-body physics [1]. This problem has re¬ 
cently attracted great interest due to experimental advances 
where interactions, mass imbalance and even the dimension¬ 
ality of the system can be controlled, thus allowing a clean 
platform to study various properties. Up to now, a handful 
of states such as the polaron [2, 3], dressed dimer or dimeron 
[4-7], and trimeron [8, 9] have been studied (for a review, see 
[10, 11]). They can be best understood as quasiparticles: few- 
body states dressed by particle-hole fluctuations of the Fermi 
sea whose vacuum properties are qualitatively preserved in the 
gas. 

Surprisingly, for the polaron, a simple variational calcula¬ 
tion [12] with only one particle-hole dressing gives very ac¬ 
curate results compared with more sophisticated Monte Carlo 
calculations [4, 13] and can be extended to study other possi¬ 
ble states of the system. For example, away from the unitary 
regime and approaching the EEC side (Fig. 1), the impurity 
binds one particle from the Fermi sea and forms a dimeron 
ground state [5, 6]. The trimeron becomes the ground state if 
the ratio of the masses of the impurity to that of the Fermi sea 
atoms is above a critical mass ratio [9] . 

The polaron, dimeron and trimeron have well-defined vac¬ 
uum limits when the density of the background gas tends to 
zero (the bare impurity, dimer and trimer) [14]. Therefore it is 
interesting to ask whether there are states which do not have a 
well-defined vacuum limit. A new such class would be bound 
states where holes play an essential role (beyond being part of 
the particle-hole dressing). For example, a bound state of the 
impurity with a hole could not be interpreted as such in the 
low density limit but its excitonic-type physics would come 
into play at nonzero densities. It is interesting to study this 
possibility in this context; excitonic states are, of course, well 
known already in charged systems, particularly in solid state 
[ 1 ]. 

In this paper, we study variationaly the simplest case of 
an impurity bound with a single hole (an exciton) within a 
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Figure 1. (color online) Schematic for equal masses in the wide 
resonance limit of the key concepts discussed in the paper. Besides 
the dimeron-hole continuum identified previously (red band spanned 
by vectors A) [15], there is also a metastable exciton-particle contin¬ 
uum (green band spanned by vectors B) and related exciton physics. 
The bottom of the A band corresponds to the dimeron energy while 
that of the B band is the exciton energy and the tops are estimated as 
being roughly the quasiparticle energies plus £p. 


two-channel model characterised by a scattering length and 
effective range in 3D and compare its energy with that of 
the polaron and dimeron (see Fig.l). We find that the exci¬ 
ton is not the ground state and calculate its decay rates to the 
polaronic and dimeronic branches, finding that it is a long- 
lived quasiparticle. An important consequence is that the exis¬ 
tence of the exciton gives rise to an exciton-particle continuum 
in inverse-rf measurements which overlaps with the already 
known dimeron-hole continuum ([11, 15]) as we see in Fig. 1. 
As a consequence, it is likely that part of the spectral weight 
of the continuum measured in that experiment and identified 
as dimeron-hole was most likely of the exciton-particle type 
instead. 
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II. MODEL 


We consider an atomic impurity of mass M immersed in a 
spin polarized Fermi gas of mass m. The interaction between 
the impurity and the background atoms is tunable through 
an s-wave Feshbach resonance described by a two-channel 
model and the interactions between the background atoms can 
be neglected at low temperatures. The Hamiltonian describing 
the system is given by (e.g., [16]) 


^=I 




- Jk + {E^ + Vq) ^k^k 


V k 


( 1 ) 


where the coupling of the impurity and a background atom 
into the closed channel molecule is described by g. Here ul^ 
(^k and ^k) creates a background atom (an impurity and a 
closed channel molecule) with wave vector k. The Wk obey the 
usual fermionic anticommutation relations {^k^^k'} = ^k' 
and the statistics of the impurity and the molecule are ir¬ 
relevant since there is at most one impurity and one closed 
channel molecule in the system. We define /2m, 

Ek = h\^/2M, £’k = /2{M^m) where Vq is the internal 

energy of the molecule. The bare parameters g and Vq of the 
Hamiltonian are related to the physical parameters, scattering 
length a and effective range of the resonance as [16] 
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where fl = mM/ (m+M) is the reduced mass. To describe the 
excitonic state of momentum |p) of the impurity in the Fermi 
sea, we use the grand canonical ensemble [ 10 ] and consider 
the variational wave function with one particle-hole dressing. 
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together with the standard ansatze for the polaron and dimeron 

[ 10 , 11 ], 


\^p)p — [ + ^0q^p+q^q + ^^q^+q-k^k^q ) 1^^) 


kq 


(4) 


\^d)p — I + ^^k^_k^k + 

V k kq 


p+q-k^k^q 


+ Y, ^kk'q^p+q_k-k'4'4^q 1^^) (5) 

kk'q / 

where |FS) is the zero-temperature noninteracting Fermi sea 
of background atoms with a Fermi energy Ep and wave vec¬ 
tor = 2m£F- Here, means that the summation of the 


hole momentum q is restricted to ^ < ^77 while the particle 
momentum is restricted to k > ^ 77 . 

The equations that determine the polaron, dimeron and ex- 
citon energies are given in the large system limit by the min¬ 
imisation of ('FjT? — £fN — L’I'F), N being the total number 
operator for background atoms. Setting the zero of the energy 
at that of the grand canonical energy of a gas of noninteracting 
atoms with the same kp, we get 
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where the various quantities in the above equations are defined 
as 


^q — ^pol ^p+q ^0 £q 

^r'q = ^pol - e^+q-k - £k + £q, (9) 


zpexc _ 77 _ e" I p^ _ p^ 

~ ^exc ^p+q ' ^q 


^qq^ — ^exc ^p+q+q' ^0 + Pq 


qq' 

776X0 _ 77 r«( 

^kqq' — ^exc “ 


£q, - £p 


and 


rexc 

qcxc _ qq^ 

^‘‘qq' “ g2 




p+q+q'-k 
1 




— fck + tq + fcq/ — £p 
1 ' 

^ v ‘^qq' ’ 


rdim _ 77 r^d i ^ 

^k ~ ^dim ~ tp-k ~ Pk~^^F 
^kq ~ ^dim “ ^p+q—k ~ ~ ^k^ ^q^ ^F 

77 dim _ T7 (yd (yU (yli I (yU 1 (y 

^kk'q - ^dim - tp+q-k-k' “ % - ^k + + £7^ 

^dim / 


f-vdim_ 

^^kq — 


"kq 


-E — 

y Z 7 dim 


1 ' 1 

7 ~ (^dim ~ ^p ~ Pf) / g ~ ^Y^ ’ 

k ^k 


k' ^kk'q 
1 

y ^ 


( 10 ) 


( 11 ) 


for polaron, exciton and dimeron respectively. Note that the 
energy equation of polaron is algebraic while that for the exci¬ 
ton and dimeron are integral equations. From the above equa¬ 
tions, we can extract the exciton, polaron and dimeron ener¬ 
gies as a function of various parameters, like, a, and M/m 
(see [5, 16]). In the next section we will consider the cases 
with p = 0. 


III. ENERGY OF THE EXCITON 

In Fig. 2 we present the energies of the exciton compared 
with that of the polaron and dimeron as a function of {kpa)~^ 
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Figure 2. (color online) The energies of the polaron, exciton and 
dimeron as a function of the interaction {kjpa)~^ with different mass 
M/m and resonance parameter kR''. Panels (al)-(a3) show the ener¬ 
gies when kpR"" = 0 at three different masses M/m = 1/5,1,5 where 
the binding energy = —h^/ma^ has been subtracted for clarity. 
The dots mark the crossing points. Panels (bl)-(b3) show the ener¬ 
gies when kpR"" = 5. 


Figure 3. Equation for the T-matrix of the scattering of an impurity 
with a hole whose resonance is the exciton. Right moving lines are 
impurity lines and left moving are background atoms (holes). The 
bold line is a polaron dressing (up to one particle-hole) and the grey 
squares are the usual impurity-background atom T-matrix in the lad¬ 
der approximation, see also Fig. (4). 


at different M/m and kpR"". Panel a(2) of Fig. 2 shows the 
calculations with equal mass {M/m =1) and in the wide reso¬ 
nance limit {kpR"" = 0), where the energies of the polaron and 
dimeron have been studied previously [5-7, 13, 16]. The en¬ 
ergies of the polaron and dimeron cross at {kpac)~^ ^ 0.85 
while that of the exciton and polaron cross at {kFac)~^ ^ 
1.3. Increasing the mass ratio M/m, as in Fig. 2 (a3) where 
M/m = 5, the critical values of {kFac)~^ for polaron-dimeron 
and exciton-polaron crossings both decrease, but the exciton 
energy is alway higher than that of the dimeron. The situation 
changes when the mass of the impurity is lighter than that of 
the atoms of the Fermi sea, as in Fig. 2 (al) where M/m =1/5 
and the exciton and dimeron energies cross at {kFac)~^ ^ 0.9, 
though the polaron energy is lower than that of the exciton 
at this scale. Fig. 2 (bl)-(b3) shows the energies of the three 
states when the resonance becomes narrow {kpR^ = 5). In this 
case, the energy curves of the dimeron and exciton become es¬ 
sentially flat around unitarity and the exciton energy is alway 
higher than that of the dimeron for ah three mass ratios of 
M/m = 1/5,1,5. From this we see that the exciton is never 
the ground state but is a resonance whose lifetime we discuss 
below. 

Because we are studying variationaly a quasiparticle which 
is not the lowest energy state (as was the case of the repulsive 
polaron in [17, 18]), we should check that ansatz (3) is in¬ 
deed describing a new quasiparticle and not some trivial com¬ 
bination of other known ones which are also included in the 
Hilbert space spanned by the ansatz. For example, it might 
describe a polaron plus a hole or (if the hrst term in (3) were 
zero) a dimeron plus two holes. The energy of the latter two 
states is approximately given by the sum of the energies of 
the quasiparticles (e.g. polaron energy plus hole energy) with 
the constraint of zero total momentum. We argue as follows: 
if we look at panels (al) and (a2) of Fig. 2, we see that, in 
the hrst case, the exciton has lower energy than the dimeron 
for values smaller than {kFa)~^ 0.9. Therefore its state for 

M/m= 1/5 and for ah kpa is not trivial combination of the 
dimeron plus two holes which is of equal or higher energy. 


In the central panel, we can similarly argue that the exciton 
is lower in energy than the state of a polaron plus a hole for 
M/m = 1 since its energy is lower than that of a polaron for 
{kFa)~^ > 1.3. It follows then that, in the {M/m^{kFa)~^) 
plane, the exciton eigenstate is not trivial combination of those 
two states everywhere by continuity. Clearly, it will also not 
be trivial combination of other states that are included in the 
Hilbert space spanned by (3) such as a “Cooper pair” (dimer 
without dressing) plus two holes or single impurity plus hole 
since these are higher in energy than the dimeron or polaron 
respectively. From this we conclude that the exciton is indeed 
a new quasiparticle. 


IV. DIAGRAMMATIC FORMULATION AND LIFETIME 
OF THE EXCITON 

The variational approach is strictly speaking valid for 
eigenstates of the system, usually the ground state. The ex¬ 
citon however is not the ground state but a resonance so, 
while the variational method can be applied in a restricted 
Hilbert space as we have argued, it might be useful to present 
its diagrammatic formulation. The resonance appears in the 
impurity-hole T-matrix at positive energy and total momen¬ 
tum p as in (3). In Fig. 3 we show the diagrams correspond¬ 
ing to ansatz (3), i.e. including at most two particle-hole lines 
at any given time. We also see that intermediate states corre¬ 
spond to the repeated scattering of a polaron and a hole. These 
diagrams represent ah possible intermediate states which span 
the two-particle hole Hilbert space of the variational ansatz 
and are therefore equivalent to it. They also illustrate that the 
exciton is a true resonance, not just a combination of a po¬ 
laron and a hole since it requires interaction between the two. 
A similar type of diagrams were considered in [19]. However 
they were interested in the equal spin population case and did 
not consider the polaronic dressing of the impurity lines. 

Since the exciton is not the ground state of the system, the 
excitonic states can be a long lived quasiparticle at best since 
they can decay to a dimeron or polaron. If the lifetime is long 
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enough, then we will be able to speak of a new quasiparticle, 
the exciton, opening a new regime of metastable physics (sim¬ 
ilarly to the repulsive polaron which is also not the ground 
state but can still be observed experimentally due to its suf¬ 
ficiently long lifetime and has actually already been studied 
experimentally by several groups recently [15, 20-22]). Thus 
the exciton is a resonance in the scattering of the impurity with 
a hole, which appears as a pole at the complex energy E -f iV in 
the scattering amplitude. The real part of the pole corresponds 
approximately to the exciton energy calculated by the varia¬ 
tional ansatz. The important question then becomes - what is 
the imaginary part of the pole, i.e., what is the lifetime of the 
exciton? As discussed above, these processes require a larger 
Hilbert space than that of the variational ansatz. 

Fig. 2 shows that all three quasiparticles can decay into 
each other. The Feynman diagrams for these processes are 
shown in Fig. 4. The decay rate F^ (decay from “a” to “j3”) 
is determined by the imaginary part of the corresponding self 
energy insertions Z (the part of the diagram without external 
legs), i.e., r = —ImZ(0,C)L)) [23]. While the decay rates be¬ 
tween polaron and dimeron have been calculated at zero [24] 
and finite momentum [25] we will focus on the decay rates of 
the exciton to the dimeron and polaron. However, for com¬ 
pleteness in Fig. 4 we also show the diagrams for the decay 
between the polaron and dimeron. 

We first calculate the decay rates from exciton to dimeron, 
in which process a zero momentum exciton decays to a 
dimeron and two holes so that the decay rate involves a sum 
over the hole final momenta. The frequency sums in the 
two diagrams can be performed by contour integration which 
amounts to evaluating the frequencies at the on-shell energies 
[24]. Using a pole expansion of the dimeron and polaron prop¬ 
agator, Z)(q, co) ZoI{co — COd — /2my) and G^(q, (o) 

Zpl{o) — (Op — 12mp), with Zp,Zd and the corre¬ 

sponding residues and effective masses, we find (from the 
imaginary part of the two self-energy insertions [23]), 


t^Exc _ 

^ Dim o 


/ 






{2ny (2;r)3 

(q + q')" 


2m*r 
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where F{q) = G'^(q, cOexc + ^q)l’Z^P and = (Oexc - tOoim- 
It is interesting to note that without the matrix element term 
[F{q)—F{q')^, the decay rate scales as ^ (Aco)^, so the 
Fermi antisymmetry when swapping the two hole momenta 
q and q' gives an additional factor (Aco)^ to the decay rate. 
This is very different from the case of decay rate between po¬ 
laron and dimeron where the Fermi antisymmetry only gives 
Am [24] due to the absence of an additional particle from 
the Fermi sea as a result of the decay. From a similar cal¬ 
culation for the decay of the dimeron to exciton, we find 

We then proceed to the calculations of the decay rates from 
exciton to polaron, in which a zero momentum exciton decays 
into a polaron, a particle and two holes (the energy threshold 
for decay into a polaron and hole is higher). The two diagrams 


in the second column of Fig. 4 yield the following expression. 


-^Exc 
- Pol 


7tg\Zp r d^q d^q' d^k 
J (2;r)3 {27tf 
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where Fi(q,q') = G^(q,C0Exc + ^g)72(q + q',fi)Exc + + 

^^,)/Zp, F2(q,q') = G''(q',fi)Exc + 4")72(q + q',<OExc + ^g + 
^^,)IZp and AcOp = (Oexc — Again, we find without the 
matrix element term the above expression would give a scal¬ 
ing r\j (Am)^-^, so the Fermi antisymmetry when swapping the 
two hole momenta q and q' also gives an additional factor 
(Am)^ to the decay rate. This may be a unique feature of the 
exciton physics in this setup. A similar analysis of the two di¬ 
agrams of the polaron to exciton decay gives 

The scalings of the decay rates among the three states are 
summarised in the last column of Fig. . These scaling laws 
with high powers are interesting theoretical results on their 
own and imply very long quasiparticle lifetimes as compared 
with the usual quadratic scaling (Am)^ of the quasiparticle 
decay rates of a normal Landau fermi liquid (for the physics 
of low power laws, such as Am or even frequency indepen¬ 
dent decay rates, see [27]), and will open a new regime of 
metastable physics in Fermi gases for experimental explo¬ 
rations [25]. Relating our scaling laws to experiments, we 
can estimate the exciton lifetime, for example at the exciton- 
polaron crossing point l/kpa ^ 1.3 of Fig. 2 (a2). At the 
crossing point, the decay of exciton to polaron will van¬ 
ish, though the exciton can decay to dimeron with an en¬ 
ergy difference of Aco/Sp 0.25. Using g 2 = —\j2nl\i^a 
[24] and assuming gi ^ g 2 , we find the prefactor of the 
scaling law by doing the integral (12) numerically, result¬ 
ing in hr rsj 2 A£f{A(o/£f)^ rsj 0.01 £77 (for a Fermi energy of 
£f = hx 31kHz, the lifetime is ^ OAms), which is the typical 
scale of the repulsive polaron as measured in [15]. Relating 
this decay rate to the corresponding energy shift of the ex- 
citonic state (iF-j ^ 0.45^77), we find hr/\E- \ ^ 0.02 ^ 1, 
which shows that the excitonic state exists as a well defined, 
metastable quasiparticle in this case. 


V. DISCUSSION 

While the excitonic state is not the ground state at least for 
some parameter values, its long lifetime adds new ingredients 
to the metastable physics previously known (e.g., repulsive 
polaron). An interesting result obtained here is that the Fermi 
antisymmetry when swapping the two hole momenta q and 
q' gives an additional factor (Am)^ to the exciton decay rates 
in contrast with the previously studied decay rates between 
polaron and dimeron where only an additional factor (Am) 
was found. 

Our study opens many interesting directions since it em¬ 
phasises hole excitations as an integral part of few-body states 
in the gas. It naturally leads to questions about states with 
more than a single hole, and to analogies with solid state sys¬ 
tems even though in the latter, the physics has qualitatively 
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Figure 4. The decay processes among the three states: polaron, dimeron and exciton, which are described by the single-particle, two-particle 
and particle-hole propagators [26]. The first (second, third) column shows the diagrams of the decay processes between exciton and dimeron 
(exciton and polaron, polaron and dimeron). The last column summarises the scaling laws for these decay processes where denotes the 
decay rate of “a” to “/3”. The thick solid (dashed, wavy) line is the polaron (exciton and dimeron) propagator while the thin solid line is the 
propagator of the Fermi background atom. • denotes the exciton-polaron coupling strength gi, ■ the polaron-dimeron coupling strength g 2 
and grey square the usual impurity-background atom T-matrix in the ladder approximation. 


different aspects due to the quasiparticle charge. An interest¬ 
ing question is whether the exciton can be the ground state in 
closely related systems. For example, in 2D [18, 28-31], inte¬ 
grals over hole momenta are more important; also, one could 
think of other states such as p-wave polarons and dimerons 
[32] or dipolar impurity systems [33] and even impurities with 
spin-orbit couplings [34] . Further studies will also be needed 


to clarify how to couple the impurity with the excitonic states 
while having negligible coupling with the dimeronic or pola- 
ronic states by inverse rf (radio-frequency) process which will 
allow us to probe the excitonic states directly. 
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